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Any of two or more two-dimensional (2D) materials with similar properties can be alloyed into a new layered 
material, namely, 2D alloy. Individual monolayer in 2D alloys are kept together by Van der Waals interactions. 
The property of multilayer alloys is a function of their layer number. Here, we studied the shear (C) and layer¬ 
breathing (LB) modes of Moq jWg 5 S 2 alloy flakes and their link to the layer number of alloy flakes. The study 
reveals that the disorder effect is absent in the C and LB modes of 2D alloys, and the monatomic chain model 
can be used to estimate the frequencies of the C and LB modes. We demonstrated how to use the C and LB 
mode frequency to identify the layer number of alloy flakes deposited on different substrates. This technique is 
independent of the substrate, stoichiometry, monolayer thickness and complex refractive index of 2D materials, 
offering a robust and substrate-free approach for layer-number identification of 2D materials. 


Monolayer (SLG) and multilayer graphenes (MLGs) have 
attracted much attention due to its extraordinary propertiesij^ 
and potentials in device applications^. Besides SLG and 
MLGs, there are a large number of other two-dimensional 
(2D) systems with interesting properties. For example, NiTe 2 
and VSe 2 are semimetals; WS 2 , WSe 2 , M 0 S 2 , MoSe 2 and 
MoTe 2 are semiconductors; h-BN, and F[fS 2 are insulators; 
NbS 2 , NbSe 2 , and TaSe 2 are superconductors; Bi 2 Se 3 and 
Bi 2 Te 3 are topological insulatorsj^ Two or more layered 
materials with similar properties can be alloyed into a new 
layered material, such as Mo;cWi_;cS 2 , Moi_j:W;cSe 2 and 
(Bii_xSb;t) 2 Te 3 alloysi^- The 2D alloys are a rich source of 
2D materials because they can exhibit tunable properties be¬ 
cause the ratio of two end compositions can be controllable^ 
The key to form a good quality of 2D monolayer is mixing the 
end compositions at atomic scale. Individual monolayers in 
the 2D alloys are kept together by Van der Waals interations. 
Similar to graphite and graphene, the property of multilayer 
alloys is a function of their layer number. How to determine 
the layer number of ultrathin 2D alloys is thus of primary im¬ 
portance for fundamental science and applications. 

Several optical techniques have been developed to iden¬ 
tify the layer number of 2D materialsNevertheless, 
most techniques rely on the specific properties of 2D mate¬ 
rials. For example, one can identify the layer number, de¬ 
noted as N, of M 0 S 2 by the intensity or peak positions of 
the cotTesponding photoluminescence (PL) peak and Raman 
modes . The mostly used technique is based upon optical 
contrast of 2D materials against dielectric substrates, such as 
a Si substrate covered with a Si02 layersj ' i^ ' - To precisely 
identify N of few-layer 2D materials, the experimental op¬ 
tical contrast must be compared with the theoretical one for 
different N. However, the theoretical calculations require pre¬ 
determined parameters, including thickness of a single layer, 
complex refractive index of the material and information of 
the substrate ,"’'^ which poses great challenge for research on 
new 2D materials in early stage. Therefore, it is desirable to 
develop a novel approach to identify the layer number of few- 
layer 2D materials, relying on as little prerequisite parameters 
as possible. 


The interlayer modes of multilayer 2D materials originates 
from the relative motions of the rigid monolayer planes them¬ 
selves, either perpendicular or parallel to their normal, such as 
the shear (C) modes and the layer breathing (LB) modes3“— 
The C and LB modes are the collective vibration modes of 
all the layers so their frequency (coc and (Olb) depends on 
N and the interlayer coupling . By probing such ultralow- 
frequency (LFLF) modes in multilayer 2D materials, it is pos¬ 
sible to identify its layer number no matter what kinds of 
substrates is used to support multilayer flakes. In common 
alloys, random distribution of atoms will result in disorder 
effect on their vibration modes in terms of peak weakening 
and broadening.— For multilayer 2D alloys exfoliated from 
the bulk layered alloys, how the alloy layer modifies their C 
and LB modes is an open issue. Whether the monatomic chain 
model (MCM) based on 2D crystals^ can be applied to the C 
and LB modes in multilayer alloys for layer-number identifi¬ 
cation is not clear. 

Here, Mo 0 . 5 W 0 . 5 S 2 alloy has been used as a template, 
to explore the possibility of applying MCM to 2D alloys. 
Moi_;(Wj:S 2 alloy will be simplified as M 0 WS 2 alloy here¬ 
after for simplicity. It is demonstrated that, by only measuring 
the C and LB modes, the layer number of M 0 WS 2 alloy flakes 
up to N=ll can be determined rapidly and nondestructively. 
Because the frequencies of the C and LB modes are linked 
with the layer number by a robust general formula based on 
MCM, this technique can be applied to ultrathin 2D materials 
on any substrate. 

Ultra-thin flakes have been mechanically exfoliated from 
bulk M 0 WS 2 alloy (purchased from 2d Semiconductors Inc.), 
and transferred onto a Si substrate covered with 90-nm Si02 
film.' Raman and PL measurements are performed at room 
temperature using a Jobin-Yvon HR800 micro-Raman system 
equipped with a liquid-nitrogen-cooled charge couple detector 
(CCD), a X 100 objective lens (numerical aperture=0.90) and 
a 1800 lines/mm grating. The excitation wavelength is 466 
nm from a Ar+ laser. At this wavelength, the 1800 lines/mm 
grating enables us to have each pixel of the charge-coupled 
detector cover 0.35 cm^'. Plasma lines have been removed 
from the laser beam, using BragGrate Bandpass Alters. Mea- 
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surements down to 5 cm^* are enabled by three BragGrate 
notch filters with optical density 3 and with FWHM=5cm^^ 
Both BragGrate bandpass and notch filters are produced by 
OptiGrate Corp. The typical laser power is about 0.4 mW, to 
avoid sample heating. 





FIG. 1: (a) Raman spectrum of bulk M 0 WS 2 alloy and bulk M 0 S 2 
crystal in the ULF and 320-450cm^* ranges, (b) ULF Raman spec¬ 
tra of bulk, IL and 2L alloy flakes under HV and VV polarization 
configurations, where the features labeled by stars are residual of the 
C mode after the spectrum subtraction, (c) PL spectra of bulk, IL 
and 2L alloy flakes. The excitation wavelength is 466nm. 


2 H-M 0 WS 2 bulk alloys are formed by stacking mono- 
layer (IL) alloy together via van der Waals interactions.- 
The IL alloy contains one MoW plane sandwiched by two 
S planes, represented as S-MoW-S. Flakes containing N S- 
MoW-S layer are denoted as NL M 0 WS 2 . Raman spectrum of 
the bulk alloy is depicted in Fig. 1(a). Since, the point group 
of the bulk alloy is Dg/,, the mode at 416cm^^ is assigned as 
as it lies between Ai^ of M 0 S 2 (409cm^^) and WS 2 
(420cm^^). There are two modes at 380 cm^* and 352 cm^^ 
assigned as MoS 2 -like and WS 2 -like respectively. The 

peaks at 347 and 365 cm^' may be the disorder-related Ra¬ 
man peaks^ The Ai^, and E^g peaks exhibit broadened pro¬ 
files due to the disorder effect in the alloy compared with that 
(gray curves) in bulk M 0 S 2 , as shown in Fig. 1(a). Only one 
ULF peak is observed at 28.8 cm^^ Considering that the LB 
mode is Raman inactive and the C mode is Raman active in 
bulk MoS 2 r"'^ this peak is attributed to the C mode in bulk al¬ 
loys. Indeed, this peak is observable under both the parallel 
(VV) and cross (HV) polarization configuration, as shown in 
Fig. 1(b), similar to the case in bulk M 0 S 2 . In contrast to 
the broadening of the high-frequency mode, full width at half 
maximum (FWHM) of the C mode in bulk alloys is about 1.1 
cm^^ and almost identical to that in bulk M 0 S 2 , as displayed 
in Fig. 1(a). It suggests that the disorder effect is absent in 


the C and LB modes of 2D alloy materials. Indeed, as the C 
mode is due to rigid-layer lattice vibrations, it is reasonable 
that the random distribution of Mo and W atoms in the MoW 
layer has little influence on the interlayer vibration modes. 

The observation of the C mode in bulk alloys allows the 
deduction of (Oc of NL flakes, thus enabling the N determina¬ 
tion of multilayer M 0 WS 2 alloys. Because the disorder ef¬ 
fect is absent for the C mode in bulk M 0 WS 2 alloys, it is 
expected that it is also absent for the C and LB modes in 
NL alloys and the MCM can be applicable for NL alloys by 
treating each S-MoW-S layer as a single balli^S Therefore,the 
symmetry of NL M 0 WS 2 alloys can be analyzed in a simi¬ 
lar way as of NL M 0 S 2 crystal.^- For the NL flake, there are 
N-1 C modes and N-1 LB modes. Based on the MCM, the 
N-dependent frequencies of both the C and LB modes exhibit 
as a fan diagrams;^ There exist a series of upper and lower 
phonon branches originating from even N. The observed up¬ 
per branches of the C and LB modes originating from 2L are 
denoted as and LB^, respectively. Similar to the case 
of NL MoS 2 <^ the N-dependent frequencies of C 2 and 
branches are, respectively, given by 

(o{C 2 ){N) = (Oc{bulk) cos {n/2N), 

(o{LB 2 ){N) = (Olb (bulk) cos {n/2N), 

where {N > 2), (Oc{bulk) and C0LB{bulk) are the C and 
LB frequency of bulk alloy, respectively. C0c{bulk) = 
V2o)c{2L), and coiBibulk) = •\/2®lb( 2L). Both a)(C^)(N) 
and (o{LB 2 )(N) stiffen with increasing Ni^ 

At first, we will demonstrate how to identify IL and 2L al¬ 
loys from the observed C0c{bulk). According to the sym¬ 
metry for 2L M 0 WS 2 , both the C and LB modes , denoted 
as Eg and Aj^, respectively. We measured the ULF modes of 
many ultrathin alloy flakes with low optical contrast. There 
exist two typical sets (FA and FB) of ultrathin alloy flakes. 
Each set exhibits an identical Raman spectral feature, as de¬ 
picted in Eig. 1(b). It is obvious that the EA flake is IL alloy 
because no any peak is present above 5cm^^ in its ULE range 
under VV and HV configurations 4^ Under the VV configu¬ 
ration, the EB flake exhibits two ULE peaks located at 20.2 
cm^^ and 35.4 cm^*, respectively, while only the peak at 20.2 
cm^^ appears under the HV configuration, and equals to ex¬ 
actly 1 /times as much as (Oc{bulk) at 28.8 cm^' as given 
byEq. Q. Thus, it is indicated that the FB flake is 2L al¬ 
loy, and the 20.2 cm^ * peak corresponds to its C mode. The 
other peak is then ascribed to LB mode with a)LB(2L)=35.4 
cm^*. Indeed, the calculated C and LB mode frequencies 
of bulk M 0 WS 2 alloys based on the MCM22 are 20.2 and 
35.6 cm^* respectively, if the rigid S-MoW-W layer takes the 
average mass of {niMo + niw)/2 -\-2ms and interlayer force 

constants^ of a|^ and in M 0 S 2 , where niMo, mw and ms 
are the masses of Mo layer in M 0 S 2 , W layer in WS 2 and S 
layer in M 0 S 2 , respectively. The intrinsic LB mode in 2L al¬ 
loy can be obtained by subtracting the spectrum under the HV 
configuration from that under the VV configuration (denoted 
as VV-HV), as shown by the dashed curve in Fig. 1(b). Based 
on Eq. ([T]i, we have COLB{bulk)=50 .1 cm^^ To further support 
the layer number identification, the photoluminescence (PL) 
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spectra of FA and FB have been measured., as demonstrated 
in Fig. 1(c). The FA flake exhibits a strong PL emission at 
663 nm and no sideband appears at its lower-energy range. 
However, the FB flake exhibits a weaker PL emission at 678 
nm and there is a low-energy sideband at 768 nm, which is 
assigned to the indirect PL emission of the 2L alloy. PL emis¬ 
sion is not observed in bulk alloys at all. All the PL features 
are analogous to those in IL, 2L and bulk MoS 2 '^ 



FIG. 2: (a) Theoretical (crosses) and experimental (diamonds) fre¬ 
quencies of C 2 and phonon branches of 2L-11L alloys. The 
three # indicate the C modes in 6L and lOL alloys with frequencies 
of 0}c{2L), as indicated by the dash-dotted lines. The dashed lines 
are guides to the eye. 

The above discussion reveals that one do apply MCM for 
the C mode in 2L alloys, thus we can apply MCM to the C 
and LB modes of NL alloys without doubt. The predicted 
frequencies of the C 2 branch based on Eq. O are plotted 
in Fig. 2(a). It shows that when N>6, the dependency of 
a)(C^)(A) on N weakens with increasing N, so the differ¬ 
ence between Co{C 2 ){N) and co{C 2 ){N — 1) are too close to 
be distinguishable. In this case, we can consider another up¬ 
per phonon branch of the C modes originating from the 6L 
alloy, which is denoted as C^. The N-dependent frequencies 
of the branch is given 

0){C^){N) = (Ocibulk) cos{3n/2N) , (2) 

where (A^ > 6). Co{C^ )(N) stiffens from (Oci^L) with increas¬ 
ing N.— The predicted frequencies of the branch are also 
summarized in Fig. 2(a) by crosses. By comparing the ex¬ 
perimental frequencies of the two branches and with 
theoretical ones, 2L-11L alloys can be easily identified. De¬ 
tailed Raman spectra related with the C modes of the identi¬ 
fied 2L-1IL alloys are depicted in Fig. 2(b). The experimental 
a)(C^)(A) (2<N<6) and (o{C'^){N) (6<N<11) are summa¬ 
rized in Fig. 2(a) by diamonds, matching theoretical results 
very well. 

Besides the C modes, there are a series of LB modes in 
NL alloys analogous to the case in NL-MoS 2 <^ Because the 
LB mode in bulk alloy is Raman inactive, the upper phonon 
branches originating from 2L and 6L cannot be observed in 
M 0 S 2 and similar 2D materials To the contrary, the lower 
phonon branches originating from 2L and 6L can be easily 


observed, which are denoted as LBj and LB^ , respectively, 
whose N-dependent frequencies are given by22 


co{LB 2){N) = C0LB{bulk)fim{nl2N), (N > 2), 

(o{LBf^){N) = (OLB{bulk)sm{37:/2N),{N >6). 

Both a)(LB 2 )(N) and a)(LBg )(N) softens from (Olb{2-L) with 
increasing Ni^ The predicted frequencies of the two branches 
are summarized in Fig. 3(a). Similar to the case of C 
modes, one can also identify N for NL alloys based on the 
N-dependent ®(LB 2 )(A) or 0}{LB^){N). Raman spectra of 
the LB modes of the identified 2L-11L alloys are depicted 
in Fig. 3(b). The experimental )(A) (2<N<11) and 

Co{LB^){N) (6<N<11) have been extracted and plotted in 
Fig. 3(a) for comparison with theoretical predictions. The 
large separation between co {LB 2 ) {N) and co {LB 2 ) (N -F 1) en¬ 
sures that one can identify N of M 0 WS 2 alloy up to 1 IL only 
by the LB 2 branch. If C 0 lb{'^L) is too small in some 2D ma¬ 
terials, (o{LB 2 ){N) may quickly approach to 5-10 cm^^ , the 
detection limit of a typical Raman system. Similar to the case 
of the C mode in M 0 WS 2 alloy, one must use the LBg branch 
to identify N of thicker 2D flakes. 



FIG. 3: (a) Theoretical (crosses) and experimental (squares) frequen¬ 
cies of LB 2 and LB^ phonon branches of 2L-1 IL alloys. The three 
# indicate the LB modes in 6L and lOL alloys with frequencies of 
C0i^{2L), as indicated by the dash-dotted lines. The features labeled 
by stars in (b) are residual signals of the C mode after spectrum sub¬ 
traction. The dashed lines are guides to the eye. 

We have demonstrated how to apply MCM for identify¬ 
ing layer-number of ultrathin alloy flakes. Although the Mo 
layer in multilayer M 0 S 2 was replaced by an alloy layer with 
random distribution of Mo and W atoms to form multilayer 
M 0 WS 2 alloy, the C and LB modes in multilayer alloys can 
be well-understood by MCM that describes the rigid interlayer 
vibrations in 2D crystals. This means that MCM is adequate 
for various 2D alloys that be combined by any two end ma¬ 
terials. In principle, the above standard identification process 
can be used for common 2D crystals. Once the C^ or LB 2 
branch is detected, layer number of ultra-thin 2D materials 
can be identified. (Oc{bulk) of most attractive 2D materials 
is larger than 20 cm^* and C0LB{bulk) is usually larger than 
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COc{bulk). In this case, layer number of ultrathin 2D materials 
up to 5L can be identified from or LBj branch with a high- 
resolution Raman system. One can further identify layer num¬ 
ber of thicker flakes by examining or LBg branch once 
they are detected. 



FIG. 4: The C modes and high-frequency modes of M 0 WS 2 alloy 
flakes on quartz, 300nm Si02/Si, and 90nm Si02/Si, which were 
measured under HV and VV configurations, respectively. The inset 
shows the corresponding optical images. The scale har in each image 
is 3/iOT. 

The C and LB vibrations are an intrinsic property of multi¬ 
layer 2D materials, which are irrelevant to any substrate. In¬ 
deed, the substrate information is not involved in the above 
discussion at all. Therefore, this layer-number identification 
method proves to be substrate-free for 2D materials. As an 
example. Fig. 4 shows the C and high-frequency modes of 
four M 0 WS 2 alloy flakes exfoliated on three types of sub¬ 
strates including quartz, 300nm Si02/Si, and 90nm Si02/Si. 


Using the aforementioned method, the layer numbers of the 
corresponding flakes have been identified to be 2L and 5L. 
Their optical contrasts in their optical images are very differ¬ 
ent from each other because of different N and substrates, ' 
as shown in the inset to Fig. 4. The high-frequency modes 
of the four alloys exhibit similar spectral features as expected. 
However, the alloy flakes with the same layer number exhibit 
the C modes with almost the same frequency, in spite of the 
significant peak intensity due to the optical interference ef¬ 
fect in the MoWS2/Si02/Si multilaye n'^''^ Alloy flakes with 
different layer number display disparate C modes, offering a 
reliable diagnostic of the layer number. Similar behavior has 
also been observed for the LB modes (not shown here). It is 
evident that the identification of the layer number is not depen¬ 
dent on the substrates. Moreover, Eqs.([T]i-© are not relevant 
to the stoichiometry, monolayer thickness and complex refrac¬ 
tive index of any 2D materials. Therefore, this method can be 
applied to other 2D materials once (Oc{bulk) or (OiB{bulk) is 
determined by Raman measurement. 

In summary, we used M 0 WS 2 as a template to study the 
N-dependent ultralow-frequency Raman spectroscopy of 2D 
alloys. The interlayer C and LB modes can be observed in 
multilayer M 0 WS 2 alloys where the alloy effect can be ig¬ 
nored. (Oc and colb in NL alloys are linked with N by ro¬ 
bust formulas based on MCM. We have demonstrated how 
to identify the layer number of NL alloy by their C or LB 
modes once (Ocibulk) or atiBibulk) is determined by Raman 
measurement. In principle, this technique can be applied to 
other 2D materials produced by micro-mechanical exfolia¬ 
tions, chemical-vapor-deposition growth or transfer processes 
on various substrates, without the need for monolayer thick¬ 
ness and complex refractive index of 2D materials and the 
supporting substrate.. 
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